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The first phase of studies of the neutrino mass and mixing is essentially over. The outcome is the discovery of 
non-zero neutrino mass and determination of the dominant structure of the lepton mixing matrix. In some sense 
this phase was very simple, and nature was very collaborative with us: Two main effects - the vacuum oscillations 
and the adiabatic conversion in matter (the MSW-effect) - provide complete interpretation of the experimental 
results. Furthermore, with the present accuracy of measurements the 3;^— mixing analysis is essentially reduced 
to the 2y consideration. I will present a concise and comprehensive description of this first phase. The topics 
include: (i) the concept of neutrino mixing in vacuum and matter; (ii) physics of the oscillations and adiabatic 
conversion; (iii) the experimental evidences of the flavor transformations and determination of the oscillation 
parameters. Some implications of the obtained results are discussed. Comments are given on the next phase of 
the field that will be much more involved. 



1. Introduction 

Recent major progress in neutrino phenomenology, 
and particle physics in general, was related to studies 
of the neutrino mass and mixing. The first phase of 
these studies is essentially over, with the main results 
being 

• discovery of non-zero neutrino mass; 

• determination of the dominant structure of the 
lepton mixing: discovery of two large mixing an- 
gles; 

• establishing strong difference of the quark and 
lepton mass spectra and mixing patterns. 

Physics of this first phase is rather simple. The 
two main effects - the vacuum oscillations [l|, Q 
and the adiabatic conversion in matter (the MSW- 
effect) [1, dj. are enough for complete interpretation 
of the experimental results. (Oscillations in matter 
appear as sub-leading statistically insignificant yet ef- 
fect for the atmospheric and solar neutrino oscillations 
in the Earth.) Furthermore, at the present level of 
experimental accuracy the three neutrino analysis is 
essentially reduced to two neutrino analyzes, and de- 
generacy of parameters is practically absent. Nature 
was very "collaborative" with us, realizing the easiest 
possibilities and disentangling an interplay of various 
phenomena. 

In a sense, we have now a "standard model of neu- 
trinos" that can be formulated in the following way: 

1) . there are only three types of light neutrinos; 

2) . their interactions are described by the Standard 
electroweak theory; 

3) . masses and mixing are generated in vacuum; 
they originate from some high energy (short range) 
physics at the electroweak or/and higher scales. 

Now the goal is to test these statements and to 
search for new physics beyond this "standard model" . 
Confirmation of the LSND result by MiniBooNE 



would be discovery of such a new physics. 

The next phase of studies will be associated to new 
generation of neutrino experiments, which will start 
in 2008 - 2010. The main objectives of this new phase 
include determination the absolute scale of neutrino 
mass and sub-dominant structures of mixing: namely, 
1-3 mixing, deviation of the 2-3 mixing from maxi- 
mal value, the CP-violation phase(s). The objectives 
include also identification of neutrino mass hierarchy 
and precision measurements of already known param- 
eters. 

The next phase will be much more involved: New 
phenomena may show up at the sub-leading level. 
More complicated formalisms for their interpretation 
are required. Complete three-neutrino context of 
study will be the must. Severe problem of degeneracy 
of parameters appears. 

In these lectures"'^ I will present a concise descrip- 
tion of the first phase of studies of neutrino masses and 
mixing. I will start by a detailed discussion of the con- 
cept of neutrino mixing in vacuum and matter. In the 
second part, the the main effects involved: the vac- 
uum oscillations, oscillations in matter and the adia- 
batic conversion, are described and physics derivation 
of all relevant formulas are given. In the third part I 
will present the experimental results and existing ev- 
idences of neutrino oscillations. For each experiment 
a simple analysis is described that allows one to eval- 
uate the neutrino parameters without sophisticated 
global fit. This consideration is aimed at convincing 



^The text presented here is partially based on lectures given 
at the Les Houches Summer School on Theoretical Physics: 
Session 84: Particle Physics Beyond the Standard Model, Les 
Houches, France, 1-26 Aug 2005, as well as on materials pre- 
pared for the TASI-06 school "Exploring New Frontiers Using 
Colliders and Neutrinos" , June 4 - 30, 2006, Boulder Colorado. 
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that indeed, we see the osciUations and and our in- 
terpretation of results in terms of the vacuum masses 
and mixing is correct. 



2. Flavors, masses and mixing 
2.1. Flavor mixing 

The flavor neutrinos, Vf = (z^e, ^'/i, f^r) are defined 
as the neutrinos that correspond to certain charge lep- 
tons: e, ^ and r. The correspondence is established 
by the weak interactions: Vi and I {I = e,^,T) form 
the charged currents or doublets of the SU2 symmetry 
group. Neutrinos, vi, V2, and v^^ with definite masses 
TOi, TO2, m3 are the eigenstates of mass matrix as well 
as the eigenstates of the total Hamiltonian in vacuum. 

The vacuum mixing means that the flavor states do 
not coincide with the mass eigenstates. The flavor 
states are combinations of the mass eigenstates: 

vi = UiiVi, l = e,fi,T, i = 1,2,3, (1) 

where the mixing parameters Uu form the PMNS mix- 
ing matrix Upmns The mixing matrix can be 
conveniently parameterized as 

UpMNS = V23i023)I~sVi3iOi3)IsVi2iOi2), (2) 

where Vij is the rotation matrix in the ij-plane, 9ij is 
the corresponding angle and Is = diag{l, 1, e*^) is the 
matrix of CP violating phase. 



2.2. Two aspects of mixing. 

A number of conceptual points can be clarified us- 
ing just 21^— mixing. Also, at the present level of ac- 
curacy of measurements the 2v— dynamics is enough 
to describe the data. For two neutrino mixing, e.g. 
i^e — fa, we have 



EE 



COS ovi + sm 0i'2 , 



cos 01^2 ~ sm av\ , 



(3) 



where Va is the non-electron neutrino state, and d is 
the vacuum mixing angle. 

There are two important aspects of mixing. The 
first aspect: according to ([3]) the flavor neutrino states 
are combinations of the mass eigenstates. Propaga- 
tion of Vf, (va) is described by a system of two wave 
packets which correspond to v\ and V2. In fig.[T^). we 
show representation of v,, and Va as the combination 
of mass states. The lengths of the boxes, cos^ Q and 
sin^ 0, give the admixtures of i^i and 1^2 in I'e and Va- 
The key point is that the flavor states are coherent 
mixtures (combinations) of the mass eigenstates. The 
relative phase or phase difference of vi and 1^2 in i'e 
as well as Va is fixed: according to ([3]) it is zero in v,. 
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a). 



b). 



Figure 1: a). Representation of the flavor neutrino states 
as the combinations of the mass eigenstates. The length 
of the box gives the admixture of (or probability to find) 
corresponding mass state in a given flavor state. (The 
sum of the lengths of the boxes is normalized to 1. b). 
Flavor composition of the mass eigenstates. The electron 
flavor is shown by red (dark) and the non-electron flavor 
by green (grey). The sizes of the red and green parts give 
the probability to flnd the electron and non-electron 
neutrino in a given mass state, c). Portraits of the 
electron and non-electron neutrinos: shown are 
representations of the electron and non-electron neutrino 
states as combinations of the eigenstates for which, in 
turn, we show the flavor composition. 



and TT in Va- Consequently, there are certain interfer- 
ence effects between vi and 1/2 which depend on the 
relative phase. 

The second aspect: the relations ([3]) can be in- 
verted: 

= cos 9 Vf, — sin 9 Va, V2 = cos 9 i'a + sin 9 . (4) 

In this form they determine the flavor composition 
of the mass states (eigenstates of the Hamiltonian), 
or shortly, the flavors of eigenstates. According to (H]) 
the probability to find the electron flavor in is given 
by cos^ 9, whereas the probability that appears as 
Va equals sin^ 9. This flavor decomposition is shown 
in fig. [TJd). by colors (different shadowing). 

Inserting the flavor decomposition of mass states 
in the representation of the flavors states, we get the 
"portraits" of the electron and non-electron neutri- 
nos fig. [TJ:). According to this figure, Ve is a system 
of two mass eigenstates that, in turn, have a com- 
posite fiavor. On the first sight the portrait has a 
paradoxical feature: there is the non-electron (muon 
and tau) fiavor in the electron neutrino! The para- 
dox has the following resolution: in the v,,- state the 
i^Q-components of vi and i>2 are equal and have oppo- 
site phases. Therefore they cancel each other and the 
electron neutrino has pure electron flavor as it should 
be. The key point is interference: the interference of 
the non-electron parts is destructive in Vg. The elec- 
tron neutrino has a "latent" non-electron component 
which can not be seen due to particular phase arrange- 
ment. However, during propagation the phase differ- 
ence changes and the cancellation disappears. This 
leads to an appearance of the non-electron compo- 
nent in propagating neutrino state which was origi- 
nally produced as the electron neutrino. This is the 
mechanism of neutrino oscillations. Similar consider- 
ation holds for the Va state. 
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2.3. Who mixes neutrinos? 

How mixed neutrino states (that is, the coherent 
mixtures on the mass eigenstates) are created? Why 
neutrinos and not charged leptons? In fact, these 
are non-trivial questions. Creation (preparation - 
in quantum mechanics terms) of the mixed neutrino 
states is a result of interplay of the charged current 
weak interactions and kinematic features of specific 
reactions. Differences of masses of the charged lep- 
tons play crucial role. 

Let us consider three neutrino species separately. 

1) . Electron neutrinos: The combination of mass 
eigenstates, which we call the electron neutrino, is pro- 
duced, e.g., in the beta decay (together with electron). 
The reason is the energy conservation: no other com- 
bination can be produced because the energy release is 
about few MeV, so that neither muon nor tau lepton 
can appear. 

2) . Muon neutrino. Almost pure v^^ state is pro- 
duced together with muons in the charged pion decay: 
7r+ — > jX^Vfi- Here the reason is "chirality suppres- 
sion" - essentially the angular momentum conserva- 
tion and V-A character of the charged current weak 
interactions. The amplitude is proportional to the 
mass of the charged lepton squared. Therefore the 
channel with the electron neutrino: 7r+ — > e'^v^ is 
suppressed as cx ml/mf^. Also coherence between i/^ 
and small admixture of i^e is lost almost immediately 
due to difference of kinematics. 

3) . Tau neutrino. Enriched i^t - flux can be ob- 
tained in the beam-dump experiments at high ener- 
gies: In the thick target all light mesons (tt, K which 
are sources of usual neutrinos) are absorbed before 
decay, and only heavy short living particles, like D 
mesons, have enough time to decay. The D mesons 
have also modes of decay with emission of z/g and 
i/^ that are chirality-suppressed in comparison with 
D tVt- Furthermore, coherence of Vf, and with 
Vr is lost due to strongly different energies and mo- 
menta. 

What about the neutral currents? Which neutrino 
state is produced in the decay in the presence of 
mixing? Z*^— interactions are flavor blind and all the 
neutrino flavors are produced with the same ampli- 
tude (rate). The only characteristic that distinguishes 
neutrinos is the mass. So, the state produced in the 
Z^-decay can be written as 

I/) = ^ W^i^i) + \^2V2) + \i>:^y:i)] (5) 

(which is also equivalent to the sum of pairs of the 
flavor states) 0]. It is straightforward to show that 
the decay rate \Z^) — > |/) is given by 

\{f\H\Z^)\^^?,\{D,v,\H\Z^)\\ (6) 

that coincides with what one obtains in the case of 
three independent decay channels. 



Do neutrinos from Z"-decay oscillate? One can 
show that oscillations can be observed in the two- 
detector experiments when both neutrinos from the 
decay are detected [6|. If a flavor of neutrino (an- 
tineutrino) is fixed, then a flavor of the accompany- 
ing antineutrino (neutrino) will oscillate with distance 
and energy. 



3. Physics effects 

3.1. To determination of oscillation 
parameters 

In the Table [J we show parameters to be deter- 
mined, sources of information for their determination 
and the main physical effects involved. In the first ap- 
proximation, when 1-3 mixing is neglected, the three 
neutrino problem splits into two neutrino problems 
and parameters of the 1-2 and 2-3 sectors can be de- 
termined independently. 

Essentially two effects are relevant for interpreta- 
tion of the present data in the lowest approximation: 

• vacuum oscillations (both averaged and non- 
averaged) 0, i2, 3]; 

• adiabatic conversion in medium 0, Q . 

A priory another effect - oscillations in matter - 
should also be used in the analysis. It is relevant for 
the solar and atmospheric neutrinos propagating in 
the matter of the Earth. It happens however, that for 
various reasons the effect is small - at {1 — 2)a level 
and can be neglected in the first approximation. 

In the case of solar neutrinos, for the preferable val- 
ues of oscillation parameters of the LMA solution (see 
below) this effect is indeed small. Furthermore, due 
to the attenuation (see below) the Earth-core effect 
is small and one can consider oscillations as ones in 
constant density. 

In the case of atmospheric neutrinos the v^— and 
Vf^~ transition probabilities driven by 1-2 mixing and 



Table I Parameters and effects. 

Parameters Source of information Main physics effects 
Ami2, Oi2 Solar neutrinos Adiabatic conversion 

and averaged vacuum 
oscillations 

KamLAND Non-averaged vacuum 

oscillations 

Am23, 023 Atmospheric neutrinos Vacuum oscillations 
K2K Vacuum oscillations 

6i3 CHOOZ Vacuum oscillations 

Atmospheric neutrinos Oscillations in matter 
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mass splitting are not small (of the order one in the 
sub-GeV range). However, due to an accidental coin- 
cidence (the fact that the ratio of the muon-to-electron 
neutrino fluxes equals 2) the effect cancels for maxi- 
mal 2-3 mixing (see below). 

Notice also that the 2v— mixing analyzes are 
enough. However, in the next order, when sub-leading 
effects are included, the problem becomes much more 
difficult and degeneracy of parameters appear. We 
will comment on this later. 



3.2. Neutrino oscillation in vacuum 

In vacuum, the neutrino mass states are the eigen- 
states of the Hamiltonian. Therefore dynamics of 
propagation has the following features: 

• Admixtures of the eigenstates (mass states) in 
a given neutrino state do not change. In other 
words, there is no vi <-> U2 transitions, vi and 1/2 
propagate independently. The admixtures are 
determined by mixing in a production point (by 
0, if pure flavor state is produced). 

• Flavors of the eigenstates do not change. They 
are also determined by 9. Therefore the picture 
of neutrino state (fig. [T]c) does not change dur- 
ing propagation. 

• Relative phase (phase difference) of the eigen- 
states monotonously increases. 

The phase is the only operating degree of freedom 
and we will consider it in details. 

Phase difference. Due to difference of masses, 
the states vi and V2 have different phase velocities 
Vphase = Ei/pi « 1 +mf/2E'^ (for ultrarelativistic 
neutrinos), so that 



Aw 



^ phase 



2E^ ' 

The phase difference changes as 

(j) = AVphaseEt. 



(7) 



(8) 



Explicitly, in the plane- wave approximation we have 
the phases of two mass states (pi = Eit — piX. Ap- 
parently, the phase difference which determines the 
interference effect one should be taken in the same 
space-time point: 



-(1)2^ AEt - Apx. 



Since p ~ \J E'^ — m^, we have 

dp . ^ dp . n 1 . ^ Am^ 
Ap = -£AE + -T-^Am^ = —AE - — - 
dE drw^ Vg Zp 



(9) 



, (10) 



where Vg = dE/dp is the group velocity. Plugging 
(fTU)) into dH) we obtain 



AEit 



2p 



(11) 



Depending on physical conditions either AE w 
or/and (t — x/vg) is small which imposes the bound on 
size of the wave packet. As a consequence, the first 
term is small and we reproduce the result ([8]). For 
stationary source one should take AE — 0. 

In general, depending on conditions of production 
and detection both quantities AE and Ap are non- 
zero. There is always certain time interval in the prob- 
lem. At, that determines (according to the uncertainty 
principle) the energy interval AE. E.g. in the case of 
solar neutrinos we know a time interval (determined 
by the time resolution of a detector) when a given neu- 
trino is detected. Furthermore, neutrino production 
processes have certain life-times, or coherence times. 
There are arguments that one should take the center 
of the wave packet where t = x/vg, or average over 
the wave packet length that leads to vanishing the 
first term in pT|) . In both cases one obtains standard 
expression for the phase. Apparently, the oscillation 
effect should disappear in the limit Am^ — 0. 

Notice that oscillations are the effect in the configu- 
ration space. The process is described by interference 
of the wave functions that correspond to the mass 
eigenstates, ipi{x,t) and ip2{x,t). Formally, we can 
perform the Fourier expansion of these wave func- 
tions considering the interference in the momentum 
representation. So, formally we can always take the 
same momenta doing then appropriate integration. 

Increase of the phase leads to the oscillations. In- 
deed, the change of phase modifies the interference: 
in particular, cancellation of the non-electron parts 
in the state produced as Ve disappears and the non- 
electron component becomes observable. The process 
is periodic: when <j) — -k, the interference of non- 
electron parts is constructive and at this point the 
probability to find Va is maximal. Later, when (p = 2tt, 
the system returns to its original state: i>{t) = v^. The 
oscillation length is the distance at which this return 
occurs: 



27r 



AttE 



AVphaseE At 



(12) 



The depth of oscillations, Ap, is determined by the 
mixing angle. It is given by maximal probability to 
observe the "wrong" flavor j/q. From the fig. Ic. one 
finds immediately (summing up the parts with the 
non-electron flavor in the amplitude) 



Ap = (2 sin 9 cos 9f = sin^ 261. 



(13) 
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Putting things together we obtain expression for the 
transition probabihty 



P = A, 



1 



cos ^— I = sin"' 20 sin^ , „ . (14) 



I. 



AE 



The oscihations are the effect of the phase increase 
which changes the interference pattern. The depth of 
oscillations is the measure of mixing. 

3.3. Paradoxes of neutrino oscillations 

A number of issues in theory of neutrino oscillations 
is still under discussion. Here I add several comments. 
Naive plane-wave description reproduces correct re- 
sult since it catches the main feature of the effect: 
phase difference change. Clearly it can not explain 
whole the picture because the oscillations are a finite 
space-time effect. 

Field theory approach provides with a consistent 
description. Oscillation experiment includes neutrino 
production in the source, propagation between the 
source and detector, detection. In this approach pro- 
duction, propagation and detection of neutrinos are 
considered as a unique process in which vi and 1/2 
are virtual particles propagating between the produc- 
tion, xp, and detection, xd, points. Propagation of 
Vi [i = 1, 2) is described by propagators Si{xD — xp). 
Notice that here there is a substantial difference from 
our standard calculations of the probabilities and 
cross-sections when we consider the asymptotic states 
and perform integration over the infinite space-time. 
The later leads to appearance of the delta- functions 
that express conservation of the energy and momen- 
tum. In the case of oscillations integration should be 
performed over finite production and detection regions 
(integration over xp and xd)- Also one should take 
into account finite accuracy of measurements of the 
energy and momenta of external particles. 

From this point of view in usual consideration we 
perform truncation oi whole process: For |a;p — a;£)| 2> 
1/Ap neutrinos can be considered as real (on-shell) 
particles with negligible corrections due to virtuality. 
Whole the process can be truncated in three parts: 
1). production; 2). propagation, as propagation of 
wave packets; 3). detection. Neutrino masses are ne- 
glected in the production and detection processes. In 
this picture, the oscillations are considered as the ef- 
fect of propagation with certain initial and final condi- 
tions that reflect process of production and detection. 
(Their effects develop over much larger space-time in- 
tervals.) Correct boundary (initial and final) condi- 
tions should be imposed. Essentially these conditions 
determine the length and shape of the wave packets. 

Let us stress again that oscillations are the finite 
space and finite time phenomenon: all the phases of 
the processes, production, propagation and detection 
occur (and should be considered) in the finite time 
intervals and finite regions of space. 



3.4. Evolution equation 

In vacuum the mass states are the eigenstates of 
Hamiltonian. So, their propagation is described by 
independent equations 



.dv., 

l—r- 

dt 



Pi 



2pi 



(15) 



where we have taken ultrarelativistic limit and omit- 
ted the spin variables that are irrelevant for the flavor 
oscillations. In the matrix form for three neutrinos 
V = {vi, V2, v^)^ , we can write 



(16) 



where Mj-^^ = diag{m\,m\,m\). Using the relation 

v = Upj^yij^gi'f we obtain the equation for the 
flavor states: 



dt 



2E 



(17) 



where = UpMNs\Mdiag\'^Upj^jj^g is the mass ma- 
trix squared in the flavor basis. In (fT7|) we have omit- 
ted the term proportional to the unit matrix which 
does not produce any phase difference and can be ab- 
sorbed in the renormalization of the neutrino wave 
functions. So, the Hamiltonian of the neutrino sys- 
tem in vacuum is 



Hn = 



|M|2 
2E 



In the 2i^ mixing case we have explicitly: 



Ar 



AE 



-cos 261 sin 26* 
sin 26 cos 20 



(18) 



(19) 



Solution of the equation (jl7[) with this Hamiltonian 
leads to the standard oscillation formula (fT4|) . 

3.5. Matter effect 

Refraction. In matter, neutrino propagation is af- 
fected by interactions. At low energies the elastic for- 
ward scattering is relevant only (inelastic interactions 
can be neglected) It can be described by the po- 
tentials Ve, Va- In usual medium difference of the 
potentials for v,, and Ua is due to the charged current 
scattering of Ve on electrons (i^ge — > Vgc) [4]: 



V ^Ve-Va^ V2G 



FUe 



(20) 



where Gp is the Fermi coupling constant and is 
the number density of electrons. The result follows 
straightforwardly from calculation of the matrix el- 
ement V = {^\Hcc\'^), where ^ is the state of 
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medium and neutrino. Equivalently, one can describe 
the effect of medium in terms of tlie refraction index: 
Href - I = V/p. 

Tlie difference of tlie potentials leads to an appear- 
ance of additional phase difference in the neutrino sys- 
tem: (jimatter = (K — Va)t. The difference of poten- 
tials (or refraction indexes) determines the refraction 
length: 



lo = 



27r 



K - K 



(21) 



Iq is the distance over which an additional "matter" 
phase equals 2it. 

In the presence of matter the Hamiltonian of system 
changes: 



V, 



(22) 



where Hq is the Hamiltonian in vacuum. Using (|T8 
we obtain (for 2j/ mixing) 



H 



m 

2E 



+ V, V^diag{V,0). 



(23) 



The evolution equation for the flavor states in mat- 
ter then becomes 



-cos 26* sin 26* 
sin 20 cos 29 



V 



vf. (24) 



The eigenstates and the eigenvalues change: 



2E' 2E 



Him, H2m- 



(25) 
(26) 



The mixing in matter is determined with respect to 
the eigenstates of the Hamiltonian in matter, and 
i'2m- Similarly to ([3]) the mixing angle in matter, 9^, 
gives the relation between the eigenstates in matter 
and the flavor states: 



COSUmVlm - 
: COs9ml'27n 



sin 9mV2m, 

- sinOmVim- 



(27) 



The angle 9m in matter is obtained by diagonalization 
of the Hamiltonian in matter 



sin^ 26'„ 



sin^ 261 



(cos 26* - 2VE/l\m'^ f + sin^ 261 



(28) 



In matter both the eigenstates and the eigenvalues, 
and consequently, the mixing angle depend on matter 
density and neutrino energy. It is this dependence 
activates new degrees of freedom of the system and 
leads to qualitatively new effects. 



Resonance. Level crossing. According to (|28p . the 
dependence of the effective mixing parameter in mat- 
ter, s\VL'29m, on density, neutrino energy as well as 
the ratio of the oscillation and refraction lengths: 



2EV 



has a resonance character. At 



(X EUe 



(29) 



Iq cos 26* (resonance condition) (30) 



the mixing becomes maximal: sin 29m =^ 1- For small 
vacuum mixing the condition (j30p reads: 

Oscillation length w Refraction length. (31) 

That is, the eigen-frequency which characterizes a sys- 
tem of mixed neutrinos, 1/1^, coincides with the eigen- 
frequency of medium, I/Iq. 

For large vacuum mixing (cos 26*12 = 0.4—0.5) there 
is a significant deviation from the equality (|3ip . Large 
vacuum mixing corresponds to the case of strongly 
coupled system for which the shift of frequencies oc- 
curs. 

The resonance condition ([50]) determines the reso- 
nance density: 



Am^ cos 20 

2E y^Gp 



(32) 



The width of resonance on the half of the height (in 
the density scale) is given by 



2Anf = 2nf tan 26*. 



(33) 



Similarly, one can introduce the resonance energy and 
the width of resonance in the energy scale. 

In medium with varying density, the layer where the 
density changes in the interval 



i^±An^ 



(34) 



is called the resonance layer. In resonance, the level 
splitting (difference of the eigenstates H2m — Him ) 
is minimal 0, Q and therefore the oscillation length 
being inversely proportional the level spitting, is max- 
imal. 



3.6. Soft neutrino masses 

One possible deviation from the standard scenario 
can be related to existence of the "soft neutrino" 
masses or situation when a part of neutrino masses 
are soft. The neutrino masses can be generated by 
some low energy physics, so that the masses change 
with energy (distance) scale; also environment effect 
on the masses becomes substantial. Recently, such a 
possibility has been considered in the context of Ma- 
VaN scenario [o^. Neutrino mass is some function of 
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a small VEV, v of some new scalar field nii, — m^{v), 
and in turn, v can depend on an environment, and 
in particular, on density of the background neutrinos 
{e.g. relic neutrinos). Another possibility is that the 
effective neutrino mass is generated by the exchange of 
light scalar boson that couples with usual matter (lep- 
tons and quarks). Scalar interactions lead to chirality- 
flip and therefore to generation ot true mass (and not 
just change of the dispersion relation as in the case of 
refraction.) Denoting the coupling constants of scalar 
boson with neutrinos and charged fermions by and 
A/ correspondingly, we find the soft mass 



(35) 



So, in the evolution equation that describes oscilla- 
tions one has 



(36) 



where ruyac is a mass generated by some short range 
physics, e.g., the electroweak scale VEV. 



3.7. Degrees of freedom 

An arbitrary neutrino state can be expressed in 
terms of the instantaneous eigenstates of the Hamil- 
tonian, vim and h'2m, as 

l^{t) = COsOal^lni + Sineal^2me"^ , (37) 



where 



m 9a = Oa(t) determines the admixtures of eigen- 
states in iy{t); 



• (pit) is the phase difference between the two 
eigenstates (phase of oscillations): 



(j){t) = I AHdt' + <j){t)i 
Jo 



(38) 



here AH = Him — H2m- The integral gives 
the adiabatic phase and (f>{t)T can be related 
to violation of adiabaticity. It may also have a 
topological contribution (Berry phase) in more 
complicated systems; 



3.8. Oscillations in matter. Resonance 
enhancement of oscillations 

In medium with constant density the mixing is con- 
stant: 9m{E,n) = const. Therefore 

• the flavors of the eigenstates do not change; 

• the admixtures of the eigenstates do not change; 
there is no ^ i^2m transitions, and h'2m 
are the eigenstates of propagation; 

• monotonous increase of the phase difference be- 
tween the eigenstates occurs: A(f)„i = {H2m ~ 

Hlm)t- 

This is similar to what happens in vacuum. The 
only operative degree of freedom is the phase. There- 
fore, as in vacuum, the evolution of neutrino has a 
character of oscillations. However, values of the oscil- 
lation parameters (length, depth) differ from those in 
vacuum. They are determined by the mixing in mat- 
ter and by the effective energy splitting in matter: 



sin^ 261 



sin^ 29,n, 



liy > I'm — 



2n 



H 



2m 



(39) 



For a given density of matter the parameters of os- 
cillations depend on the neutrino energy which leads 
to characteristic modification of the energy spectra. 
Suppose a source produces the i>e- flux Fq(E). The 
flux crosses a layer of length, L, with a constant den- 
sity He and then detector measures the electron com- 
ponent of the flux at the exit from the layer, F{E). In 
fig-Hlwe show dependence of the ratio F[E)/Fq{E) on 
energy for thick and thin layers. The oscillatory curve 
is inscribed in to the resonance curve (1 — sin^ 29m) ■ 
The frequency of the oscillations increases with the 
length L. At the resonance energy, the oscillations 
proceed with maximal depths. Oscillations are en- 
hanced in the resonance range: 

E = Er± AEr, AEr = Er tan 261 = E% sin 26*, 

(40) 

where E^}^ = Am^ /2^/2GFne. Notice that for E > 
Er, matter suppresses the oscillation depth; for small 
mixing the resonance layer is narrow, and the oscilla- 
tion length in the resonance is large. With increase of 
the vacuum mixing: Er and AEr — > E^. 

The oscillations in medium with nearly constant 
density are realized for neutrinos of different origins 
crossing the mantle of the Earth. 



• 9m{ne{t)) determines the flavor content of the 
eigenstates: (veWim) ^ cos 9m, etc.. 

Different processes are associated with these three 
different degrees of freedom. 



3.9. MSW: adiabatic conversion 

In non-uniform medium, density changes on the way 
of neutrinos: Ue = n^it). Correspondingly, the Hamil- 
tonian of system depends on time, H = H{t), and 
therefore, 
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Figure 2: Resonance enhancement of oscillations in 
matter with constant density. Shown is a dependence of 
the ratio of the final and original fluxes, F/Fo, on energy 
{x = Iv/lo oc E) for a thin layer, L — lo/n (left panel) and 
thick layer L — IO/q/ti" (right panel), lo is the refraction 
length. The vacuum mixing equals sin^ 20 = 0.824. 



(i) the mixing angle changes during propagation: 

= d,n{ne{t)); 

(ii) the (instantaneous) eigenstates of the Hamilto- 
nian, vim and J^2m, are no more the "eigenstates" of 
propagation: the transitions vim ^ t^2m occur. 

However, if the density changes slowly enough the 
transitions vi„i V2m can be neglected. This is 
the essence of the adiabatic condition: and V2m 
propagate independently, as in vacuum or uniform 
medium. 

Evolution equation for the eigenstates. Adiabaticity. 
Let us consider the adiabaticity condition. If external 
conditions (density) change slowly, the system (mixed 
neutrinos) has time to adjust this change. 

To formulate this condition let us consider the evo- 
lution equation for the eigenstate of the Hamiltonian 
in matter. Inserting J^/ = U{9m)vm in to equation for 
the flavor states (IMt we obtain 



.dvm 
I — — 
dt 




As follows from this equation for the neutrino eigen- 
states [1, [13], \0m\ determines the energy of transition 
Vim ^ vira and |7?2m " ii\m\ givBS the energy gap 
between levels. 
If [3 



il2m — H 



l7n 



« 1, 



(42) 



the off-diagonal terms can be neglected and system 
of equations for the eigenstates split. The condition 
(j42p means that the transitions <-> z^2m can be ne- 
glected and the eigenstates propagate independently 
(the angle 9a (|57)) is constant). 

For small mixing angles the adiabaticity condition 
is crucial in the resonance layer where (i) the level 
splitting is small and (ii) the mixing angle changes 
rapidly. If the vacuum mixing is small, the adiabatic- 
ity is the most critical in the resonance point. It takes 
the form Q 



Arfl > Ir, 



(43) 



where In = Z^/ sin 20 is the oscillation length in 
resonance, and Arji ~ urJ {dn^/ dr)Rt&Ti29 is the 
spatial width of resonance layer. 

MSW- effect. Dynamical features of the adiabatic 
evolution can be summarized in the following way: 

• The flavors of the eigenstates change according 
to density change. The flavor composition of the 
eigenstates is determined by 9m{t). 

• The admixtures of the eigenstates in a propagat- 
ing neutrino state do not change (adiabaticity: 
no Uim V2m transitions). The admixtures are 
given by the mixing in the production point, 9^. 

• The phase difference increases; the phase veloc- 
ity is determined by the level splitting (which in 
turn, changes with density (time)). 

Now two degrees of freedom become operative: the 
relative phase and the flavors of neutrino eigenstates. 
The MSW effect is driven by the change of flavors 
of the neutrino eigenstates in matter with varying 
density. The change of phase produces the oscillation 
effect on the top of the adiabatic conversion. 

Let us derive the adiabatic formula H, [ll|, [T^ . 
Suppose in the initial moment the state Vf, is produced 
in matter with density tiq. Then the neutrino state 
can be written in terms of the eigenstates in matter 
as 



cos 9^\l^ir 



sin 0° 



V2r. 



(44) 



(41) 



where 9^ = 9mino) is the mixing angle in matter in 
the production point. Suppose this state propagates 
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non-oscillatory conversion 



n° > n 



•Vl' 



conversion + oscillations 




oscillations + small matter effect 



Figure 3: Adiabatic evolution of neutrino state for three 
different initial condition (rig). Shown are the neutrino 
states in different moments of propagation in medium 
with varying (decreasing) density. The yellow vertical 
line indicates position of resonance. The initial state is Vc, 
in all the cases. The sizes of the boxes do not change, 
whereas the flavors (colors) follow the density change. 



adiabatically to the region with zero density (as it 
happens in the case of the Sun). Then, the adiabatic 
evolution will consists of transitions vim vi , V2m — > 
z^2, and no transition between the eigenstates occurs, 
so the admixtures are conserved. As a result the final 
state is 

|z.(t)) ^cosCl'^i) +smCe*^<*V2>, (45) 

where (p is the adiabatic phase. The survival proba- 
bility is then given by 



P=\{yeHt))\ 



(46) 



Plugging \v{t)) (|45|) and given by ([3]) into this 
expression and performing averaging over the phase 
which means that the contributions from and \ v2) 
add incoherently, we obtain 

P = (cos6lcos6l,"„)2 (sin6'sin6l°„)2 

= sin^ 61 + cos 261 cos^ dl, . (47) 

This formula gives description of the solar neutrino 
conversion with accuracy 10"^, that is, corrections 
due to the adiabaticity violation are extremely small 

Physical picture of the adiabatic conversion. Ac- 
cording to the dynamical conditions, the admixtures 
of eigenstates are determined by the mixing in neu- 
trino production point. This mixing in turn, depends 
on the density in the initial point, n^, as compared 
to the resonance density. Consequently, a picture of 
the conversion depends on how far from the resonance 
layer (in the density scale) a neutrino is produced. 

Three possibilities relevant for solar neutrino con- 
version are shown in fig. [3l The state produced as 
propagates from large density region to zero density. 
Due to adiabaticity the sizes of boxes which corre- 
spond to the neutrino eigenstates do not change. 



1) . TT-e S> n^: production far above the resonance 
(the upper panel). The initial mixing is strongly sup- 
pressed, and consequently, the neutrino state, v,,, con- 
sists mainly of one {v2m) eigenstate, and furthermore, 
one flavor dominates in this eigenstate. In the reso- 
nance (its position is marked by the yellow line) the 
mixing is maximal: both flavors are present equally. 
Since the admixture of the second eigenstate is very 
small, oscillations (interference effects) are strongly 
suppressed. So, here we deal with the non-oscillatory 
flavor transition when the flavor of whole state (which 
nearly coincides with V2m) follows the density change. 
At zero density we have V2m = ^^2, and therefore 
the probability to find the electron neutrino (survival 
probability) equals 

P = \{yeHtW ~ \{^e\^2,n{t))? - \{l^e\l^2)? ~ sinH . 

(48) 

This result corresponds to 9^ = 7r/2 in formula (|47)) . 

The value of final probability, sin^ 6, is the feature 
of the non-oscillatory transition. Deviation from this 
value indicates a presence of oscillations. 

2) . rig > nf^: production above the resonance (mid- 
dle panel). The initial mixing is not suppressed. Al- 
though i>2rn is the main component, the second eigen- 
state, I'lm, has an appreciable admixture; also the fla- 
vor mixing in the neutrino eigenstates is significant. 
So, the interference effect is not suppressed. As a re- 
sult, here an interplay of the adiabatic conversion and 
oscillations occurs. 

3) . rig < n^: production below the resonance 
(lower panel). There is no crossing of the resonance 
region. In this case the matter effect gives only cor- 
rections to the vacuum oscillation picture. 

The resonance density is inversely proportional to 
the neutrino energy: oc 1/E. So, for the same 
density profile, the condition 1) is realized for high 
energies, the condition 2) is realized for interme- 
diate energies and condition 3) - for low energies. 
As we will see all three case realize for solar neutrinos. 

The adiabatic transformations show universality: 
The averaged probability and the depth of oscilla- 
tions in a given moment of propagation are deter- 
mined by the density in a given point and by initial 
condition (initial density and fiavor). They do not de- 
pend on density distribution between the initial and 
final points. In contrast, the phase of oscillations is 
an integral effect of previous evolution and it depends 
on a density distribution. 

The universal character of the adiabatic conversion 
can be further generalized in terms of variable [5| 



y = 



_ nf 



(49) 



which is the distance (in the density scale) from the 
resonance density in the units of the width of reso- 
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nance layer psp . In terms of n the conversion pattern 
depends only on the initial value y^. 



p 

" ^ 












O.EV 




0.4 ' 




0.2 





Figure 4: The dependence of the average probability 
(dashed hne) and the depth of oscillations (p"''^, p"^" 
solid lines) on y for yo — —5. The resonance layer 
corresponds to y = 0. For tan'^ d = 0.4 (large mixing 
MSW solution) the evolution stops at j// = 0.47. 

In fig. m we show dependences of the average prob- 
ability, P, and depth of oscillations determined by 
pniax pmm^ y rj,^^ probability itself is the 
oscillatory function which is inscribed into the band 
shown by the solid lines. The average probability is 
shown by the dashed line. The curves are determined 
by the initial value yo only. In particular, there is no 
explicit dependence on the vacuum mixing angle. The 
resonance is at y and the resonance layer is given 
by the interval y — — 1 1. The figure corresponds 
to yo = —5, i.e., to production above the resonance 
layer; the oscillation depth is relatively small. With 
further decrease of yo, the oscillation band becomes 
narrower approaching the line of non-oscillatory con- 
version. For zero final density we have 



1 



Vf 



tan 26* 



(50) 



So, the vacuum mixing enters the final condition. For 
the best fit LMA point, y/ = 0.45 - 0.50, and the 
evolution should stop rather close top the resonance. 
The smaller mixing the larger final yf and the stronger 
transition. 



3.10. Adiabaticity violation 

In the adiabatic regime the probability of transition 
between the eigenstates is exponentially suppressed 
Pi2 ^ and 7 is given in (|42l) [n','lZ|. One 

can consider such a transition as penetration through 
a barrier of the height — Him by a system with 
the kinetic energy dOm/dt. 

If density changes rapidly, so that the condition 
is not satisfied, the transitions vim V2m become ef- 
ficient. Therefore the admixtures of the eigenstates 



in a given propagating state change. In our pictorial 
representation (fig. 3) the sizes of boxes change. Now 
all three degrees of freedom of the system become op- 
erative. 

Typically, adiabaticity breaking leads to weakening 
of the fiavor transition. The non-adiabatic transitions 
can be realized inside supernovas for the very small 
1-3 mixing. 



4. Determination of the oscillation 
parameters 

4.1. Solar neutrinos 

Data. Data analysis is based on results from 
the Homestake experiment [l3 |. Kamiokande and Su- 
perKamiokande Il5| . from radiochemical Gallium ex- 
periments SAGE [161, Gallex [l^ and GNO flS^ and 
from SNO [19] . The information we have collected can 
be described in three-dimensional space: 

1. Type of events: ve scattering (SK, SNO), CC- 
events (CI, Ga, SNO) and NO events (SNO). 

2. Energy of events: radiochemical experiments in- 
tegrate effect over the energy from the threshold to 
the maximal energy in the spectrum. Also NC events 
are integrated over energies. The CC events in SNO 
and ve events at SuperKamiokande give information 
about the energy spectrmn of original neutrinos. 

3. Time dependence of rates (searches for time 
variation of the fiux). 

Evidence of conversion. There are three types of 
observations which testify for the neutrino conversion: 

1 ) . Deficit of signal which implies the deficit of the 
electron neutrino flux. It can be described by the ratio 
R EE iV°''VAf^^*^, where Af'S'SM jg ^he signal predicted 
according to the Standard solar model fluxes [201 ■ The 
deficit has been found in all (but SNO neutral current) 
experiments. 

2) . Energy spectrum distortion - dependence of the 
suppression factor on energy. Indirect evidence is pro- 
vided by comparison of the deficits in experiments sen- 
sitive to different energy intervals: 

Low energies (Ga) : R = 0.5 — 0.6 (51) 
High energies {CI, SK, SNO) : R « 0.3. (52) 

So, the deficit increases with neutrino energy. 

3) . Smallncss of ratio of signals due to charged cur- 
rents and neutral currents [191 1: 



CC 

— = 0.340 ± 0.023 (stat. 



tVofii^y^t.) . (53) 



The latter is considered as the direct evidence of the 
fiavor conversion since NC events are not affected by 
this conversion, whereas the number CC events is sup- 
pressed. 
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All this testifies for the LMA MSW solution. 
Till now there is no statistically significant observa- 
tions of other signatures of the conversion, namely, 

- distortion of the boron neutrino spectrum: up 
turn at low energies in SK and SNO (significant ef- 
fect should be seen below 5-7 MeV); 

- day-night effect (recall that SK agrees with pre- 
dictions however significance is about Icr); 

- semiannual time variations on the top of an- 
nual variations (due to eccentricity of the Earth orbit). 

Physics of conversion [21] . Physics can be described 
in terms of three effects: 

1) . Adiabatic conversion (inside the Sun); 

2) . Loss of coherence of the neutrino state (on the 
way to the Earth); 

3) . Oscillations of the neutrino mass states in the 
matter of the Earth. 

According to LMA, inside the Sun the initially pro- 
duced electron neutrinos undergo the highly adiabatic 
conversion: —> cos^JJ^i^i + sm9^i'2, where 9^ is the 
mixing angle in the production point. On the way 
from the central parts of the Sun the coherence of 
neutrino state is lost after several hundreds oscillation 



lengths j2lj, and incoherent fluxes of the mass states 
i>i and 1/2 arrive at the surface of the Earth. In the 
matter of the Earth vi and 1/2 oscillate partially regen- 
erating the z^e-flux. With regeneration effects included 
the averaged survival probability can be written as 

P = sin^ 9 + cos^ 6'5"2° cos 26*12 - cos 29'^^frea- (54) 

Here the first term corresponds to the non-oscillatory 
transition (dominates at the high energies), the second 
term is the contribution from the averaged oscillations 
which increases with decrease of energy, and the third 
term is the regeneration effect, with the regeneration 
factor, freg defined as 



fr. 



P,. 



(55) 



Here P2e is the probability of 1/2 — >■ i^e transition in 
the matter of the Earth (without oscillations in mat- 
ter: P2e — sm^9). At low energies P reduces to the 
vacuum oscillation probability with very small matter 
corrections. 

There are three energy ranges with different features 
of transition: 

1. In the high energy part of spectrum, £^ > 10 
MeV (x = 1„/Iq > 2), the adiabatic conversion with 
small oscillation effect occurs. At the exit, the result- 
ing averaged probability is slightly larger than sin^ 9 
expected from the non-oscillatory transition. With 
decrease of energy the initial density approaches the 
resonance density, and the depths of oscillations in- 
creases. 

2. Intermediate energy range E ^ {2 — 10) MeV 
{x — 0.3 — 2 ) the oscillation effect is significant. The 
interplay of the oscillations and conversion takes place. 



For E ^ 2 MeV neutrinos are produced in reso- 
nance. 

3. At low energies: E <2 MeV (x < 0.3), the vac- 
uum oscillations with small matter corrections occur. 
The averaged survival probability P « 0.5 sin^ 29 
is given by approximately the vacuum oscillation 
formula. 

Inside the Earth. Entering the Earth the state z/2 
(which dominates at high energies) splits in two mat- 
ter eigenstates: 



1^2 



(56) 



It oscillates regenerating partly the z/g-flux. In the 
approximation of constant density profile the regener- 
ation factor equals 

7 EV 
f = 0.5-^ sin^ 261 = - — - sin^ 26*. (57) 
Iq Am"^ 

Notice that the oscillations of 1/2 are pure matter effect 
and for the presently favored value of Am^ this effect 
is small. According to ([ST]) , freg oc l/Am^ and the 
expected day-night asymmetry of the charged current 
signal equals 



AdN = freg/P ^ (3 - 5)% 



(58) 



Apparently the Earth density profile is not con- 
stant and it consists of several layers with slow density 
change and jumps of density on the borders between 
layers. It happens that for solar neutrinos one can 
get simple analytical result for oscillation probability 
for realistic density profile. Indeed, the solar neutrino 
oscillations occur in the so called low energy regime 
when 



2EV(x) 

e= . V « 1' 
Am^ 



(59) 



which means that the potential energy is much smaller 
than the kinetic energy. For the LMA oscillation pa- 
rameters and the solar neutrinos: e{x) = (1 — 3)-10^^. 
In this case one can use small parameter e{x) ([59)1 to 
develop the perturbation theory (2^ . The following 
expression for the regeneration factor, freg, has been 
obtained [H [H 



dxV (x) sin (fimix X f) . (60) 



freg = - sin^ 261 



Here xq and x / are the initial and final points of prop- 
agation correspondingly, and (j)rn{x — > x/) is the adi- 
abatic phase acquired between a given point of tra- 
jectory, x, and final point, Xf. The latter feature has 
important consequence leading to the attenuation ef- 
fect - weak sensitivity to the remote structures of the 
density profile when non-zero energy resolution of de- 
tector is taken into account. On the other hand freg 
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can be strongly affected by some relatively small struc- 
tures near the surface of the Earth. 

Another insight into phenomena can be obtained 
using the adiabatic perturbation theory which leads 
to 



Jreg — 



2E sin-^ 29 



f E AT/, sin f. (61) 

j=0...n-l 



Here and are the phases acquired along whole 
trajectory and on the part of the trajectory inside the 
borders j. This formula corresponds to symmetric 
profile with respect to the center of trajectory. Using 
(|6T|) one can easily infer the attenuation effect. The 
formula reproduces precisely the results of exact 
numerical calculations. Notice that the adiabatic 
perturbation theory is relevant here because the adi- 
abaticity is fulfilled within the layers and maximally 
broken at the borders. 

Determination of the solar oscillation parameters. 
Knowledge of the energy dependence of the adiabatic 
conversion allows one to connect the oscillation pa- 
rameters with observables immediately. 

1). Determination of the mixing angle. To explain 
stronger deficit at higher energies one needs to have 
6 < 7r/4 or sin^ 6* < 1/2. Furthermore, using the fact 
that Ph > sin^ 9 and Pi < 0.5 sin^ 29 we find 



> 



sm 



1 



Pi ~ 0.5sin2 26' 2cos2 6l' 



(62) 



where on the RHS we have taken the asymptotic val- 
ues of the survival probability. Consequently, 



sm 



< 1 - 



2Ph 



0.1 ±0.2. 



(63) 



The ratio of CC to NC events determines the sur- 
vival probability: 



P = s\n^ 9 + cos29{cos^ ft" 



CC 



(64) 



For high energies and without Earth matter regener- 
ation effect P = sin^ 9. Since no significant distor- 
tion of the energy spectrum is seen at SK and SNO 
the Boron neutrino spectrum should be in the flat 
part (bottom of the "suppression pit"). In this re- 
gion the deviation from asymptotic value is weak. For 
Am^ ~ 8 • 10~^ eV^ the averaged oscillation effect is 
about 10%. Therefore 



CC 

sin2 « 0.9— « 0.31 



(65) 



2). Determination of Am^. Value of suppression 
in the Gallium experiments, Pi, implies that the pp- 
spectrum is in the vacuum dominated region, whereas 
stronger suppression of SK and SNO signals (together 



with an absence of distortion) means that the boron 
neutrino flux is in the matter dominated region. So, 
the transition region should be Etr ^ (1 — 4) MeV. 
On the other hand the expression for the middle en- 
ergy of the transition region equals (it corresponds to 
neutrino production in resonance) 



Am^ cos 29 



(66) 



where Vprod is typical potential in the neutrino pro- 
duction region in the Sim. From (j66p we obtain 



Am^ = 



2EtrVr 



tr ^ prod 



COS 26* 



(67) 



which gives l\m^ = (3 — 15) • 10^^ eV^ in the correct 
range. 

Another way to measure Am^ is to study the 
high energy effects: according to LMA the splitting 
Am^ is restricted from below by the increasing 
day-night asymmetry and from above by absence of 
the significant up turn of spectrum at low energies. 

New SNO results are expected from the third 
(last) phase of the experiment that employs the '^He- 
counters for neutrons. The counters provide with a 
better identification of the NC-events and therefore 
preciser measurements of the CC/NC ratio, and the 
9i2 angle (combination of cos"* ^13 sin^ 9i2 in the three 
neutrino context). BOREXINO should start measure- 
ments soon [25} . 

The SAGE calibration result is about 2a below 
expectation [2^. That may testify for lower cross- 
section and therefore higher pp— flux at the earth due 
to larger survival probability. That produces some 
tension in the fit of the solar data [l^]. Another pos- 
sibility proposed recently is that the reduced calibra- 
tion result is due to short range oscillations to sterile 
neutrinos [28|. 

Searches for time variations and possible periodicity 
in the solar neutrino data are continued 12411 . 



4.2. KamLAND 

KamLAND (Kamioka Large Anti-neutrino detec- 
tor) is the reactor long baseline experiment [2§| . Few 
relevant details: Ikton liquid scintillator detector sit- 
uated in the Kamioka laboratory detects the antineu- 
trinos from surrounding atomic reactors (about 53) 
with the effective distance (150 - 210) km. The clas- 
sical reaction of the inverse beta decay, Vf,p — > e^n, is 
used. The data include 

(i) the total rate of events; 

(ii) the energy spectrum (fig. [5]); 

(iii) the time dependence of the signal which is due 
to variations of the reactors power. (Establishing the 
correlation between the neutrino signal and power of 
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reactors is important check of the whole experiment). 
In fact, this change also influences the oscillation effect 
since the effective distance from the reactors changes 
{e.g., when power of the closest reactor decreases). 
In the oscillation analysis the energy threshold E > 



2.6 MeV prompt 
analysis threshold 



KainLAND data 
best-fit oscillation 
best-fit decay 
best-fit decoherence 




40 50 60 

L^/E- (km/MeV) 

Figure 5: The L/E distribution of events in the 
KamLAND experiment; from [2^ . 

The physics process is essentially the vacuum oscil- 
lations of De. The matter effect, about 1%, is negligi- 
ble at the present level of accuracy. 
The evidences of the oscillations are 
1). The deficit of the number of the De events 



N, 



obs 



258 



N, 



expect 



365.2 ±23.7 



0.7 



(68) 



for E > 2.6 MeV. 

2). The distortion of the energy spectrum or L/E 
dependence (when some reactors switch off the ef- 
fective distance changes). Notice that an absence of 
strong spectrum distortion excludes large part of the 
oscillation parameter space. 

Oscillation parameters are related to the observ- 
ables in the following way. The main features of the 
L/E dependence are maximum at {L/E)max — 32 
km/MeV (phase </> — 2tt) and 

16, 48 km/MeV (0 = 7r,37r). They fit well the ex- 
pected oscillation pattern. Taking the first maximum 
we find 



An 



{L/E)r 



8 • 10^^ eV^ 



The deficit of the signal determines (for a given 
Am^) the value of mixing angle: 



sin-" 2e^ 



l-Ru 



(69) 



where the averaged over the energy interval oscil- 
latory factor can be evaluated for the KamLAND 
detector as (sin^ (j)) ^ 0.6. Notice that sensitivity to 
mixing angle is not high at present. 



Extracted values of the oscillation parameters are 
in a very good agreement with those obtained from 
the solar neutrino analysis. This comparison implies 
the CPT conservation. 

Combined analysis of the solar neutrino data and 
the KamLAND can be performed in assumption of 
the CPT conservation. The mixing angle is mainly 
determined by the solar neutrino data, whereas Am^ 
is fixed by the KamLAND. New complete calibration 
of the detector will allow to improve sensitivity to 1-2 
mixing (30j . 

Comparison of results from the solar neutrinos and 
KamLAND open important possibility to check the 
theory of neutrino oscillation and conversion, test 
CPT, search for new neutrino interactions and new 
neutrino states. 



4.3. Atmospheric neutrinos 

Experimental results. The atmospheric neutrino 
flux is produced in interactions of the high energy cos- 
mic rays (protons, nuclei) with nuclei of atmosphere. 
The interactions occur at heights (10 - 20) km. At 
low energies the fiux is formed in the chain of decays: 
TT — > /ii'^, /i — !■ eVf,Vy^. So, each chain produces 2v^i and 
\ve, and correspondingly, the ratio of fluxes equals 



2. 



(70) 



With increase of energy the ratio increases since the 
lifetime acquires the Lorentz boost and muons have 
no time to decay before collisions: they are absorbed 
or loose the energy. As a consequence, the fiux of the 
electron neutrinos decreases. 

In spite of the long term efforts, still the predicted 
atmospheric neutrino fluxes have large uncertainties 
(about 20% in overall normalization and about 5% in 
the so called "tilt" parameter which describes the un- 
certainty in the energy-dependence of the flux). The 
origin of uncertainties is twofold: original flux of the 
cosmic rays and cross sections of interactions. 

The recent analyzes include the data from Baksan 
telescope, SuperKamiokande [HI, [H, MACRO [H], 
SOUDAN [341 . The data can be presented in the three 
dimensional space which includes 

- type of events detected: e-like events (showers), 
/i-like events, multi-ring events, NC events (with de- 
tection of 7r°), r— enriched sample of events. 

- energy of events: widely spread classification in- 
cludes the sub-GeV and multi-GeV events, stopping 
muons, through-going muons, etc.. 

- zenith angle (upward going, down going, etc). 
Now MINOS experiment [sl] provides some early 

information on oscillation effects for the atmospheric 
neutrinos and antineutrinos separately. 
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The evidences of the atmospheric neutrino oscilla- 
tions include: 

1). Smallness of the double ratio of numbers of ji- 
like to e-like events ISlf : 



R 



'IK 



Nf" /N, 



th ■ 



(71) 



The ratio weakly depends on energy: it slightly in- 
creases from sub-GeV to multi-GeV range (as ex- 
pected): 

i?^/e = 0.658 ±0.016 ±0.035 (subGeV) 

R^,,e = 0.702t°;°^2 ± 0.101 (multiOeV + PC). (72) 

Apparently in the absence of oscillations (or other 
non-standard neutrino processes) the double ratio 
should be 1. The smallness of the ratio testifies for 
disappearance of the j/^— flux. 

2). Distortion of the zenith angle dependence of the 
ji -like events (see fig. E]). Global characteristic of this 
distortion is the up-down asymmetry defined as 



A 



up/ down 



Ndo 



(73) 



Due to complete up-down symmetric configuration for 
the production, in the absence of oscillations or other 
non-standard effects the asymmetry should be absent: 



A 



up /do 



1. 



The zenith angle dependence for different types of 
events in different ranges of energies is shown in fig. [S] 
from [3l[. The zenith angle of the neutrino trajectory 
is related to the baseline L as L = D cos O^- So, study- 
ing the zenith angle distributions we study essentially 
the distance dependence of the oscillation probability. 

Substantial distortion of the zenith angle distribu- 
tion is found. The deficit of numbers of events in- 
creases with decrease of cosdz and reaches about 1/2 
in the upgoing vertical direction for multi-GeV events. 
The distortion increases with energy. Correspond- 
ingly, the up-down asymmetry increases with energy: 
In contrast to the /i-like, the e-like events distri- 
bution does not show any anomaly. Though one can 
mark some excess (about 15%) of the e-like events in 
the sub-GeV range (upper-left panel of fig. 0. 

3) . Appearance of the r-like events (2.4(t effect) 

4) . The i/i?— dependence shows the first oscilla- 
tion minimum (fig.[7|)- 

In the first approximation all these data can be con- 
sistently described in terms of the i/^ — Vr vacuum os- 
cillations. Notice that for pure 2i/ oscillations of this 
type no matter effect is expected: the matter poten- 
tials of j/^ and Vr are equal. In the context of three 
neutrino mixing, for non-zero values of sin ^13 the mat- 
ter effect should be taken into account for the — 
channel. 



multi-ring 
Sub-GeV |i-lii<e 
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Figure 6: The zenith angle distribution of the 
atmospheric fi-\ike events in different energy ranges; from 
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Figure 7: L/E distribution of the atmospheric fi like 
events; from [3^ . The solid line corresponds to the 
oscillation fit. 



As we marked above, the probabilities of Vg and 
oscillations in matter of the Earth driven by the "so- 
lar" parameters Am^^ and sin^ 2^12 are large and even 
matter-enhanced in the sub-GeV range. However, ob- 
servable effects of these oscillations are suppressed by 
factor 

(rcos^ 023 - 1), (74) 



where the ratio r is defined in eq. ([70|). In the sub- 
GeV range r « 2 and for maximal 2-3 mixing effects 



IPM-LHP06-19May 



IPM School on Lepton and Hadron Physics, Tehran, Iran, May 15-20, 2006 



cancel. With increase of neutrino energy r increases, 
however the probabihties are suppressed by matter 
effect. 

So, in the first approximation a unique description 
in terms of — Vt oscillation is valid for different 
types of events and in a very wide range of energies: 
from 0.1 to more than 100 GeV. 

Determination of the atmospheric neutrino oscilla- 
tion parameters. Let us describe how the oscillation 
parameters can be immediately related to observables. 
We will use here the interpretation of the results in 
terms of 2z^-oscillations v^^ — v^-. 

The most clean way to determine parameters is to 
use the zenith angle distribution of the multi-GeV [i- 
like events. As follows from fig. [S] for the down-going 
muons, co^Oz ~ 0.5-^1, the oscillation effects are neg- 
ligible (good agreement with the no-oscillation predic- 
tions). For the up-going muons, cos 6*2 ^ —0.5 — 1, 
there is already the averaging oscillation effect. Tran- 
sition region corresponds to the horizontal events with 
cosOz ~ 0.0 0.2. For these events the baseline 
L — 500 km should be comparable with the oscilla- 
tion length L K li,, so that 

AnEjnulU-GeV 



(75) 



Taking {1-2) GeV we find Am^ = (1 - 4) • 10"^ 
eV^. An uncertainty in the neutrino direction and 
the fact that distance strongly depends on cos 9z in 
the horizontal direction lead to the uncertainty in the 
determination of the atmospheric Am^. 

For the upward-going /i-like events the oscillations 
are averaged (no dependence of the suppression factor 
on cosOz), so that N°'"'{up)/N*''{up) = 1 - sin^ 26*. 
This allows us to determine the mixing angle: 

sin^ 26* = 2 [1 - N°'"' (up) /TV*'' (up)] . (76) 

From the fig. [6] TV"'"* (up) /iV*'' (up) w 0.5, and conse- 
quently, sin^ 261 = 1 {N^^{up) « N°^'{up)). 

Other independent determinations are possible: in 
the sub-GeV range the zenith angle dependence is 
weak because of strong averaging effect: (i) the oscil- 
lation length is shorter and therefore the oscillations 
develop already for large part of the downgoing events; 
(ii) the angle between neutrino and detected muon is 
very large and directionality is essentially lost. So, 
taking the deficit of the total number of events we 
obtain 



sin^ 20 > 2[1 - R], 



(77) 



where equality corresponds to the developed oscilla- 
tions for all directions. From fig. [5] we find R = 0.67, 
and therefore sin^ 29 > 0.7. 

To determine mixing angle one can use also the dou- 
ble ratio. As follows from ([TTjl 



sin^ 261 ^ 



1-R 



where (sin^ 4'l'^)z is the averaged over the energy and 
zenith angle oscillatory factor. For multi-GeV events 
(sin^ (l)/2), = 0.20 - 0.25 and therefore from ^ we 
obtain sin^ 29 1. 

The most precise determination of Am^ follows 
from the L/E - dependence of the events (fig. [7]) which 
is considered as the direct observation of the neutrino 
oscillations - oscillatory effect 36|. In the first oscil- 
lation minimum - dip in the survival probability the 
phase of oscillations equals = tt. Therefore 



Am^ = 



27r 



(79) 



/i/e 



;sin2 0/2),' 



(78) 



{L/E)dip 

From fig. m {L/E)dip ^ 500 km/GeV. This gives 
immediately Am^ = 2.5-10"^ eV^. 

Results of the 3:/— analysis from [sj and [111 show 
one important systematic effect: the shift of 2-3 mix- 
ing from maximal one when the effect of the 1-2 sector 
is included. Also whole allowed region is shifted. Even 
larg er deviation of sin^ 023 from 0.5 has been found in 
p38| . Essentially this result is related to the excess of 
the e-like events in the sub-GeV range. 

4.4. K2K 

The Vp— beam with typical energies E — (0.5 — 3) 
GeV was created at KEK and directed to Kamioka. 
Its interations were detected by SuperKamiokande 
[soj . The baseline is about 250 km. The oscillations of 
muon neutrinos, v^, have been studied by com- 

parison of the detected number and the energy spec- 
trum of the /x-like events with the predicted ones. The 
predictions have been made by extrapolating the re- 
sults from the "front" detector to the Kamioka place. 
The front detector similar to SK (but of smaller scale) 
was at about 1 km distance from the source. 

The evidences of oscillations were (i) the deficit of 
the total number of events: 107 events have been ob- 
served whereas 15ltio have been expected; (ii) the 
spectrum distortion (fig. [8]). Searches for the j/^ — > Ve 
oscillations gave negative result. 

The data are interpreted as the non-averaged vac- 
uum oscillations v^ — Vt. 

The energy distribution of the detected /z— like 
events show an evidence of the first oscillation dip 
at £; ~ 0.5 GeV (see fig. ^. This allows one 
to evaluate Am^. Using the relation (I79p with 
L/E = 250km/0.5GeV = 500 km/GeV (apparently 
the same as in the atmospheric neutrino case), we ob- 
tain Am? — 2.5 • 10^'^ eV^ in perfect agreement with 
the atmospheric neutrino result. (In fact the data 
stronger exclude other values of Am^ than favor the 
best one.) 

The substantial oscillation suppression is present 
in the low energy part of the spectrum {E < 1 
GeV) only. Therefore the deficit of events ~ 0.67 
corresponds to large or nearly maximal mixing. 
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Figure 8: The energy spectrum of events in the K2K 
experiment, from [33|. Also shown are equally 
normalized fit curves with oscillations (solid) and without 
oscillations (dotted) 
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Figure 9: The energy spectrum of events in the MINOS 
experiment, from [401 ]. Shown also are the expected 
spectrum without oscillations and the best fit to 
experimental result. 



4.5. MINOS 

MINOS (Main Injector Neutrino Oscillation 
Search) is the long baseline experiment "from Fer- 
milab to SOUDAN". NuMI (Main Injector) beam 
consists, mainly, of v^'s with energies (1 - 30) GeV 
and the flux-maximum at ~ 3 GeV. There are two 
detectors - steel-scintillator tracking calorimeters. 
The near detector is at the distance 1 km from 
the injector with mass 1 kton and the far detector 
(SOUDAN mine) has the baseline 735 km and 
mass 5.4 kt. The first result corresponds to exposure 
1.27- 10^" protons on target. 215 neutrino events have 
been observed below 30 GeV, whereas 336.0 ± 14.4 
events were predicted on the basis of measurements 
in the near detector [40| . 

Evidence of oscillations consists of (i) deficit of the 
detected number of events - disappearance of the v^- 
flux, and (ii) distortion of energy spectrum fig. [51 
The relative suppression increases with decrease of en- 
ergy (although there is large spread of points), and the 
strongest suppression is in the bins (1 - 2) GeV. 

The dominant effect is the non-averaged vacuum 
l^^Jl, — Vr oscillations. The m atter effect is negligible. 
Taking E' ~ 1.5 GeV as the energy of the first oscilla- 
tion dip (minimum) we find Am^ ~ (2.5 — 3.0) • 10^^ 
eV^. Suppression in this bin is consistent with maxi- 
mal mixing. Using the total deficit of events one can 
evaluate the mixing angle more precisely. 

Essentially we observe the high energy part of the 
first oscillation deep starting from minimum which is 
consistent with maximal suppression. This is enough 
to make rather precise determination of Am^. De- 
tailed statistical fit gives Am^s = 2.74+"-2g x 10^^ 
eV2 (68 % C.L.) and sin^ 26*23 > 0.87 (68 % C.L.) 
with the best fit value sin^ 26*23 = 1- This is in a very 
good agreement with the atmospheric neutrino and 
K2K results. 

Comment. Simple relations we have presented in 



sect. 4.1 - 4.5 allow us to understand where sensitivity 
to different parameters comes from. These relations 
are embedded in precise statistical analysis. They al- 
low us to control the outcome of this analysis, under- 
stand uncertainties and give confidence in the results 
of more sophisticated analysis. 

They show robustness of the results and their inter- 
pretation. 

4.6. 1-3 mixing: effects and bounds 

The direct bounds on 1-3 mixing are obtained in 
the CHOOZ experiment [4l|. This is the experiment 
with a single reactor, single detector and the baseline 
about 1 km. The expected effect is the vacuum non- 
average oscillations with survival probability given by 
the standard oscillation formula: 

P= 1 -sin^ 26*13 sin^^. (80) 

The baseline is comparable with the half-oscillation 
length: For the best fit value of Am^ from the atmo- 
spheric neutrino studies and E ^ 2 MeV the oscilla- 
tion length equals ~ 2 km. 

The signature of the oscillations consists of distor- 
tion of the energy spectrum described by ([50)1 . No 
distortion has been found within the error bars which 
put the limit 

sin^ 6*13 < 0.04, (90% C.L.) (81) 

for Am^i = 2.6 • 10^"^ cV^. In the atmospheric neu- 
trinos the non-zero 1-3 mixing will lead to oscillations 
of the electron neutrinos. One of the effects would be 
i^fi t'e oscillations in the matter of the Earth. The 
resonance enhancement of oscillations in neutrino or 
antineutrino channels should be observable depending 
on the type of mass hierarchy. That can produce an 
excess of the e-like events mostly in multi-GeV range 
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where the mixing can be matter enhanced. No sub- 
stantial effect is found. Notice that in the analysis [3^ 
- the best fit value sin is non-zero due to some dis- 
tortion of the zenith angle dependence, in the multi- 
GeV range. 

In solar neutrinos, the non-zero 1-3 mixing leads 
to the averaged vacuum oscillations with small oscil- 
lation depth. The effect is reduced to change of the 
overall normalization of the flux. The combined analy- 
sis of all solar neutrino data leads to zero best-fit value 
of 1-3 mixing. The CC/NC ratio at SNO and Gallium 
results (which depend on the astrophysical uncertain- 
ties less) give sin^ 6*13 = 0.017 ± 0.026 (that indicates 
a level of sensitivity of existing observations). 

In contrast, ^13 can produce leading effects for su- 
pernova electron (anti) neutrinos. 



4.7. Degeneracy of oscillation 
parameters and global fits 

In the previous sections we have analyzed various 
data in the 21/ context. Essentially the 3j^ system 
splits in to two sectors: "solar" sector probed by the 
solar neutrino experiments and KamLAND, and the 
"atmospheric" sector probed by the atmospheric neu- 
trino experiments K2K and MINOS. This is justified 
if the 1-3 mixing is zero or small and if in the atmo- 
spheric sector studies the effect of 1-2 sector can be 
neglected. That could happen, e.g., because in the 
specific experiments the baselines are small or the en- 
ergies are large, so that oscillation effects due to 1-2 
mixing and 1-2 split have no time (space) to develop. 

In the next phase of studies when sub-leading ef- 
fects, e.g., induced by sin^ia, become important the 
splitting of problem into two sectors is not possible. 
At this sub-leading level the problem of determination 
of the neutrino parameters becomes much more com- 
plicated. 

In the tabic M wc indicate relevant parameters for 
different studies. The same observables depend on 
several parameters, so that the problem of degener- 
acy of the parameters appears. In such a situation 
one needs to perform the global fit of all available 



Table II Experiments and relevant oscillation 
parameters. 



oscillation parameter constraints 



Experiments 



Parameters of Parameters of 
leading effects sub-leading effects 
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Solar neutrinos, Am 
KamLAND 

Atmospheric neutrinos Am23, O23. 
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CHOOZ Amis, 

MINOS Ami 
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Figure 10: The results of global analysis for 1-2 and 
2-3 mass splits and mixings; from |38l] . 



data. The advantages are (1) no information is lost; 
(2) dependence of different observables on the same 
parameters is taken into account; (3) correlation of pa- 
rameters and their degeneracy is adequately treated. 

There are however some disadvantages. In partic- 
ular, for some parameters the global fit may not be 
the most sensitive method, and certain subset of the 
data can restrict a given parameter much better {e.g., 
ATO23 in atmospheric neutrinos). 

In fig. [To] we show the results of the global fit of os- 
cillation data performed in [sll before MINOS results. 
MINOS shifts the allowed region and the best fit point 
of ATO23 to larger values. With earlier MINOS result 



Am^3 



2.6 • 10~3 eV"' is foimd in 43] as the best 
fit value. 

Results of global fits of the other groups (see 
113) nil) agree very well. Different types of experi- 
ments confirm each other: KamLAND confirms so- 
lar neutrino results, K2K - the atmospheric neutrino 
results etc.. Furthermore, unique interpretation of 
whole bulk of the data in terms of vacuum masses 
and mixings provides with the overall confirmation of 
the picture So, the determination of the parameters 
is rather robust, and it is rather non-plausible that 
future measurements will lead to significant change. 

The most probable values of parameters equal 

Am?2 = (7.9 - 8.0) • 10"^ eV^ (82) 

sin^eiiz = 0.310- 0.315, (83) 

Am^a = (2.5 - 2.6) • 10"^ eV^ (84) 

sin^ 623. = 0.44 - 0.50. (85) 



Slightly smaller value of 1-2 mixing, sin^ 612 = 0.30, 
has been found in 14311. 
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The parameter which describes the deviation of the 
23 mixing from maximal equals 

D23 = 0.5 - sin^ 6*23 = 0.03 - 0.06. (86) 

For 1-3 mixing we have 

sin^ 023 = 0.00-0.01, (1(7 = 0.011-0.013). (87) 

The ratio of mass squared differences important for 
theoretical implications equals 



= 0.031 - 0.033. 



5. Neutrino mass and flavor spectrum 
5.1. Spectrum 

Information obtained from the oscillation experi- 
ments allows us to make significant progress in recon- 
struction of the neutrino mass and flavor spectrum 
(Fig. [11]). 



For charge leptons the corresponding ratio is 0.06, and 
even stronger hierarchies are observed in the quark 
sector. 

2. There is the bi-large or large-maximal mixing 
between the neighboring families (1 - 2) and (2 - 3). 
Still rather significant deviation of the 2-3 mixing from 
the maximal one is possible. 

3. Mixing between remote (1-3) families is weak. 
To a good approximation the mixing matrix has the 

so-called tri-bimaximal form [44 




2 V2 
-1 V2 V3 
1 -V2 Vs 



(90) 

where is the maximal (7r/4) rotation in the 2-3 
plane and sin^ ^12 — 1/3. Alternatively, the mixing 
can be expressed in terms of the quark-lepton com- 
plementarity (QLC) relations [45| : 

"lepton mixing = bi — maximal mixing — CKM". 

(91) 

Possible realizations are 



UpMNS - 
UpMNS 



UbrnU, 



UcKMUbm 



(92) 



where Uckm is the quark mixing matrix and Ubm is 
the bi-maximal mixing matrix: 



bra 



1^23 "^12 




(93) 



Figure 11: Neutrino mass and flavor spectra for the 
normal (left) and inverted (right) mass hierarchies. The 
distribution of flavors (colored parts of boxes) in the 
mass eigenstates corresponds to the best-fit values of 
mixing parameters and sin^ 6113 — 0.05. 

The unknowns are: 

(i) admixture of Vf. in 1/3, Uez', 

(ii) type of mass spectrum: hierarchical, non- 
hierarchical with certain ordering, degenerate, which 
is related to the value of the absolute mass scale, mi] 

(iii) type of mass hierarchy (ordering): normal, in- 
verted (partially degenerate); 

(iv) CP-violation phase 5. 



Both the tri-bimaximal mixing and the QLC-mixing 
agree with the experimental data within Icr (46j . 

5.2. Absolute scale of neutrino mass 

Direct kinematic methods. Measurements of the 
Curie plot of the decay at the end point - give 
me < 2.05 eV (95%) (Troitsk) after "anomaly" sub- 
traction [13, and me < 2.3 eV (95%), (Mainz) [H]. 
Future KATRIN experiment [i^ aims at one order of 
magnitude better upper bound: rUe < 0.2 eV (90%). 
The discovery potential is estimated so that the 
positive result — 0.35 eV can be established at 5a 
(statistical) level. 



Information described in the previous sections can 
be summarized in the following way. 

1. The observed ratio of the mass squared differ- 
ences (|88p implies that there is no strong hierarchy of 
neutrino masses: 



m2 
m3 



> 



' Am^2 



0.18 ±0.02. 



(89) 



'23 



From oscillation experiments we get the lower 
bound on mass of the heaviest neutrino: 



mil > 



Am2,„=0.04eV (95%). 



(94) 



In the case of normal mass hierarchy mh — and 
in the inverted hierarchy case m/i = mi w m2. 
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Neutrinoless double beta decay. The rate neutrino- 
less double beta decay is determined by the effective 
Majorana mass of electron neutrino 



i(j>{k) 



(95) 



r oc rrigg. Here (j){k) is the phase of the k eigenvalue. 

The best present bound on is given by the 
Heidelberg-Moscow experiment: ruee < 0.35 — 0.50 
eV 50] . Part of collaboration claims an evidence of a 
positive signal [5l|, [s^l ■ The Heidelberg-Moscow col- 
laboration searched for the mode of the decay 



76 



Ge 



.76 



Se 



+ e + 



with the end point Qee — 2039 keV. The total statis- 
tics collected from 5 enriched Ge— detectors is 71.7 kg 
yr. The peak at the end point of the total energy- 
spectrum of two electrons has been found and inter- 
preted in [ill as due to the neutrinoless double beta 
decay. Number of events in the peak gives the half- 
lifetime 



ri/2 = 1.19-lO^V, 3ct range : (0.69-4.18) -10^' y. 

(97) 

The significance of the peak depends on model of 
background and quoted by the authors as 4.2cr. There 
is a number of arguments pro and contra of such in- 
terpretation. 

If the exchange of light Majorana neutrino is the 
dominant mechanism of the decay, the measured life- 
time corresponds to the effective mass of the Majorana 
neutrino 



^25 



nie 



0.44 eV, 3cr range : (0.24 - 0.58) eV. (98) 



The H-M positive result would correspond to strongly 
degenerate neutrino mass spectrum. That, in turn, 
implies new symmetry in the leptonic sector. 

Other groups do not see signal of the 00oi, decay 
though their sensitivity is somehow lowe r Issl . Is^ . Issf . 
New experiment with '^^Ge, GERDA f56|, will be 
able to confirm the H-M claim in the first phase, and 
in the case of negative result, strongly restrict it in 
future measurements. 

Cosmological bound. Analysis of the cosmological 
data that includes CMB observations, SDSS of galax- 
ies, Lyman alpha forest observations and weak lensing 
lead to the upper bound on the sum of neutrino masses 



^m, < 0.42 eV (95% C.L.) 



(99) 



|57| (see also [5al) which corresponds to mo < 0.13 
eV in the case of a degenerate spectrum. An even 
stronger bound, ELi ™^ < 0-17 eV (95% C.L.) |5g| 



was established after publication of WMAP3 results. 
This limit disfavors a strongly degenerate mass spec- 
trum, and the positive claim of observation of neutri- 
noless double beta decay. Combining the cosmological 
and oscillation bounds, we conclude that at least one 
neutrino mass should be in the interval 



m - (0.04 - 0.10) eV (95% C.L.). 



(100) 



In future, the weak lensing will allow to perform di- 
rect measurements of clustering of all matter and not 
just luminous one. This will improve the sensitivity 
down to ^ - mi ~ 0.03 eV. 



(96) 5.3. To the new phase of the field 



In what follows we summarize the parameters, 
physics goals and physics reach of the next generation 
(already approved) LBL experiments. In each case 
we give the baseline, L, the mean energy of neutrino 
{El,) and the goals. All the estimations are given for 
the 90% C.L.. 

1) . T2K ("Tokai to Kamioka"): JPARK -> Su- 
perKamiokande |63|. This is the accelerator off-axis 
experiment on searches for — + and v^^ — > v^, 
oscillations; parameters of the experiment: L = 295 
km, {Ey) — 0.7 GeV. The goal is to reach sensitivity 
to the i^e~ appearance which will allow to put the 
bound sin^ ^13 < 0.005 (or discover the 1-3 mixing 
if the angle is larger), to measure 2-3 mass split and 
mixing with accuracies 6{/S.m\.^) ^ 0.1 meV, and 
^(sin^ 26*23) = 0.01 near the maximal mixing. The 
latter corresponds to (5(sin^ ^23) — D23 — 0.05. If 
1-3 mixing is near the present bound the hope is to 
get some information about the mass hierarchy. The 
measurements will start in 2009. 

2) . NOiyA: Fermilab Ash River [6l|. This is 
also the accelerator off-axis experiment on — > 
and —f V), oscillation searches. Parameters: 
L = 810 km, (E^) = 2.2 GeV. The objectives include 
the bound on 1-3 mixing sin^ ^13 < 0.006, precise 
measurement of Arn^^, and possibly, determination 
of the mass hierarchy. Start: 2008 - 2009. 

3) . Double CHOOZ reactor experiment (g^l will 
search for 9^ —^ oscillation disappearance. Two 
detectors setup will be employed. Parameters: 
L = 1.05 km (far detector), {E^,) = 0.004 GeV, 
L/E = 250 km/GeV; the goal is to put the bound 
sin2e'i3 < 0.005- 0.008. Start: 2008; results: 2011. 

4) . Daya Bay [Hf reactor experiment will search 
for — !■ Vf, oscillation disappearance with multi- 
detector setup: Two near detectors and one far 
detector with the baseline 1600 - 1900 m from reactor 
cores are proposed. The goal is to reach sensitivity 
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s\T? 013 < 0.0025 or better. Start: 2010. 

To large extend results from these experiments will 
determine further (experimental) developments. 

5.4. Expecting the supernova neutrino 
burst 

Detection of the Galactic supernova can substan- 
tially contribute to determination of the neutrino 
parameters and reconstruction of the neutrino mass 
spectrum. In particular this study will contribute to 
determination of the 1-3 mixing and type of the neu- 
trino mass hierarchy. 

In supernovas one expects new elements of the fla- 
vor conversion dynamics. Whole level crossing 
scheme can be probed and the effects of both MSW 
resonances (due to Am^2 and Am^g) should show up. 
Various effects associated to the 1-3 mixing can be 
realized, depending on value of ^13. The SN neutri- 
nos are sensitive to sin^ 6*13 as small as 10~^. Studies 
of the SN neutrinos will also give an information the 
type of mass hierarchy [H, H^, IH, [13] ■ The small 
mixing MSW conversion can be realized due to the 
1-3 mixing and the "atmospheric" mass split Am^g. 
The non-oscillatory adiabatic conversion is expected 
for sin^ > 10~^. Adiabaticity violation occurs if 
the 1-3 mixing is small sin^ ^13 < 10^"^. 

Collective flavor transformation effects due to the 
neutrino-neutrino scattering (flavor exchange phe- 
nomenon) can be important in the central parts (out- 
cide neutrino spheres) of the collapsing stars (gH] . 

Another possible interesting effect is related to 
shock wave propagation. The shock wave can reach 
the region of the neutrino conversion, p ~ 10* g/cc, 
after tg = (3 — 5) s from the bounce (beginning of the 
v— burst) (69j . Changing suddenly the density pro- 
file and therefore breaking the adiabaticity, the shock 
wave front influences the conversion in the resonance 
characterized by Am^g and sin^ 6'i3, ifsin2 6'i3 > 10"^ 

Monitoring the shock wave with neutrinos can shed 
some light on the mechanism of explosion (6^ . [toI . FtH . 
Ill- 

5.5. LSND result and new neutrinos 

LSND (Large Scintillator Neutrino Detector) col- 
laboration studied interactions of neutrinos from Los 
Alamos Meson Physics Facility produced in the decay 
chain: tt"*" — > + ly^, fJ^'^ — *■ e"*" -|- i/g + i^/j- The ex- 
cess of the (e+ -I- n) events has been observed in the 
detector which could be due to inverse beta decay: 
I'e+P ^ + n [7^. In turn, could appear due to 
oscillations —^ in the original beam. 

Interpretation of the excess in terms of the Vp_ — Pe 
oscillations would correspond to the transition proba- 



bility 

P = (2.64 ± 0.76 ±0.45) • 10"^. (101) 

The allowed region is restricted from below by Am^ > 
0.2 eV2. 

This result is clearly beyond the "standard" 3v pic- 
ture. It implies new sector and new symmetries of the 
theory. 

The situation with this ultimate neutrino anomaly 
[73I I is really dramatic: all suggested physical (not re- 
lated to the LSND methods) solutions are strongly 
or very strongly disfavored now. At the same time, 
being confirmed, the oscillation interpretation of the 
LSND result may change our understanding the neu- 
trino (and in general fermion) masses. 

Even very exotic possibilities are disfavored. An 
analysis performed by the KARMEN collaboration 
[t^ I has further disfavored a scenario [t^ in which the 
Pe appearance is explained by the anomalous muon 
decay Ve^iC^^ {i = e, /x, r}. 

The CPT- violation scheme [7y| with different mass 
spectra of neutrinos and antineutrinos is disfavored by 
the atmospheric neutrino data [t^ . No compatibility 
of LSND and "all but LSND" data have been found 
below 3ct [73. 

The main problem of the (3-1-1) scheme with 
Ato^ ~ 1 eV^ is that the predicted LSND signal, 
which is consistent with the results of other short base- 
fine experiments (BUGEY, CHOOZ, CDHS, CCFR, 
KARMEN) as well as the atmospheric neutrino data, 
is too small: the i^^ — > probability is about 3(T 
below the LSND measurement. 

Introduction of the second sterile neutrino with 
Am^ > 8 eV^ may help It has been shown ^ 
that a new neutrino with Am^ ~ 22 eV^ and mix- 
ings C/e5 = 0.06, = 0.24 can enhance the pre- 
dicted LSND signal by (60-70)%. The (3 -I- 2) scheme 
has, however, problems with cosmology and astro- 
physics. The combination of the two described so- 
lutions, namely the 3-1-1 scheme with CPT- violation 
has been considered [Slj. Some recent proposals in- 
cluding the mass varying neutrinos MaVaN Q and 
decay of heavy sterile neutrinos [s^l also have certain 
problems. 

MiniBooNE [s^ is expected to clarify substantially 
interpretation of the LSND result. The MiniBooNE 
searches for v,, appearance in the 12 m diameter tank 
filled in by the 450 t of mineral oil scintillator and 
covered by 1280 PMT. The flux of muon neutrinos 
with the average energy (iJ^) « 800 MeV is formed in 
TT decays (50m decay pipe) which are in turn produced 
by 8 GeV protons from the Fermilab Booster. The 541 
m baseline is about half of the oscillation length for 
Am^ ~ 2 eV^. In 2006 the experint operates in the 
antineutrino channel — > 

Of course, confirmation of the LSND result (in 
terms of oscillations) would be the most decisive 
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Figure 12: The region of oscillation parameters selected 
by LSND result versus sensitivity of the MiniBooNE 
experiment; from [s^ . 



though the problem with background should be scru- 
tinized. The negative result still may leave an am- 
biguous situation. In fig. [T^lthe sensitivity limits and 
discovery potential of MiniBooNE are shown. 



6. Conclusion 

After the first phase of studies of neutrino mass and 
mixing we have rather consistent picture: interpreta- 
tion of all the results (except for LSND) in terms of 
vacuum mixing of three massive neutrinos. Two main 
effects (consequences of mixing) are important for the 
interpretation of results at the present level of accu- 
racy: the vacuum oscillations and the adiabatic con- 
version in matter (the MSW-effect). The oscillations 
in matter give sub-leading contributions, at (1 — 2)a 
level, to the solar and atmospheric neutrino observ- 
ables. 

There are unknown yet parameters and their deter- 
mination composes a program of future phenomeno- 
logical and experimental studies. Next phase of the 
field, study of sub-leading effects, will be much more 
involved. 

The main theoretical challenge is to understand 
what is behind the observed pattern of neutrino 
masses and mixing (as well as masses and mixings 
of other fermions). What is the underlying physics? 
Clearly there is a strong difference of the quark 
and lepton mixing patterns. The data hint the tri- 
bimaximal scheme of mixing with possible implica- 
tions of new "neutrino symmetries" , or alternatively 
to the quark-lepton complementarity that hints cer- 
tain quark-lepton symmetry and unification. Are the 



tri-bimaximal or QLC relations real (follow from cer- 
tain principles) or simply accidental? 

It may happen that something important in prin- 
ciples and context is still missed. The key question 
is how far we can go in this understanding using 
our usual notions of the field theory (or the effective 
field theory) and in terms of symmetries, various 
mechanisms of symmetry breaking, etc.! The hope is 
that neutrinos will uncover something simple and illu- 
minating before we will be lost in the string landscape. 
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